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Hartree-Fock and coupled-cluster calculations have been 
performed for cubic AgCl and for AuCl having a cubic or the 
observed structure with space group lAi/amd. Cohesive en- 
ergies and lattice constants are in excellent agreement with 
experiment for AgCl; for AuCl we find good agreement, and 
the experimental structure is correctly predicted to be lower 
in energy than the cubic one. Electron-correlation effects on 
lattice constants are very large, of up to 0.8 A for cubic AuCl. 
We especially discuss the strength of the closed-shell interac- 
tions, and for the first time a quantitative analysis of the 
so-called "aurophilic" Au(I)-Au(I) interaction is presented in 
solids. 



I. INTRODUCTION 

In recent years, the role of attractive closed-shell in- 
teractions in. inorganic chemistry has been intensively 
investigateda and the-.tt;sults have been summarized in 
several review articlesH'tl. Both structural, spectroscopic, 
and energetic evidence exist. The structures exhibit short 
secondary M-M distances (M = Cu(I), Ag(I), Au(I), 
T1(I), etc.) and inwards bending of the primary L-M-L' 
bonds. The spectroscopic evidence incorporates Raman 
frequencies of the M-M stretching mode. The interac- 
tion energies can be directly measured as activation ener- 
gies by temperature-dependent NMR and from optically 
monitored chemical equilibria. 

Ab-initio theoretical studies exist for intramolecular in- 
teractions and for dimers. They invariably yield repulsive 
interactions at Hartree-Fock level and reproduce the at- 
traction at correlated leyel. Relativistic effects enhance 
the interaction for golcu. It should|-|bc added that the 
free (Au^)2 dimcr remains repulsivccl. In this sense the 
ligands are essential. 

In this article, we investigate these interactions for the 
first time in a solid, using quantum-chemical ab-initio 
methods (Hartree-Fock in combination with the coupled- 
cluster approach). We consider the systems AgCl and 
AuCl. While the former has rocksalt structure, the latter 
one exhibits a rather conuplicated tetragonal structure, 
with space group lAi/amaB, which is thcrmody^namically 
stable in a region between 90° C and 342° CU. This is 



just an example of a system with relatively short Au-Au 
distances. For the sake of comparison, we performed cal- 
culations for both the cubic and the experimental struc- 
ture, in the case of AuCl. 



II. THE METHOD 

The calculations consist of two parts. First, Hartree- 
Fock (HF) calculations on thCpSolid are performed with 
the program package CRYSTALETtEl. After that, electron 
correlation for the extended systems is included by means 
of a Bethe-Goldstone-like expansion, in terms of correla- 
tion contributions frcpj localized orbital groups ('imethod 
of local increments" Eil). In a series of studiesliil^EJ, it 
has been shown that this method is capable of yielding 
results of high quality. For a detailed account of the 
method as applied to ionic compounds, see Refs. ^jl4[ 
(In Ref. [l|, we also showed how to determine van der 
Waals Cg parameters with this method.) It should be 
mentioned that we determined individual increments in 
finite model systems (embedded clusters), using an ionic 
embedding, with cationic pseudopotentials in the first 
sphere and point charges beyond. Such an embedding 
was applied to all systems considered in the present pa- 
per, including AuCl in the experimental structure (which 
is not fully ionic, cf. below). In order to check the va- 
lidity of our description, we compared the increments 
taken from different clusters (with two or three explic- 
itly described ions) and various representations of the 
Madelung field. The experimental structure turned out 
to be more sensitive with respect to the number of point 
charges than the rocksalt structure, so that a larger num- 
ber is necessary in the former case. However, the results 
with 8x8x8 unit cells were nearly identical (up to a 
few /LtH) to those of calculations where only 6x6x6 
unit cells were represented by point charges, which en- 
sures convergence of the results with the larger set. A 
more severe problem was the transferability of the incre- 
ments which was less good in the observed AuCl structure 
than in the cubic one. As correlation scheme, we chose 
the (size-consistent) coupled-cluster approach with sin- 
gle and double substitutions (CCSD), and the CCSD(T) 
scheme including triples by perturbation theory. 

For Ag and Au, we used energy-consistent relativis- 
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tic 19-valence-electron pseudopotentialsE-^, in combina- 
tion with the corresponding optimized {8s7p6d) valence 
basis sets (uncontracted) and augmented with 3/ and 2g 
functions (Ag //g-exponents are 2.54, 0.73, 0.20 / 1.58, 
0.50; Au //5-exponents are 1.41, 0.47, 0.15 / 1.20, 0.40). 
Results for the neutral atoms and singly charged cations, 
with these pseudopotcntials and basis sets, are shown 
in Table |; the ionization potentials turn out to be in 
good agreement with experiment. For CI, an all-electron 
description was chosen; the basis set used for the cor- 
relation calculatioiis, is a [6s5p3d2f] augmented valence 
triple zeta basis setlla already employed in earlier calcula- 
tions on alkali halidesEil. Finally, single- valence-electron 
pseudopotentialstZI for Ag and Au were used as part of 
the embedding, in our calculations, in order to model 
the cage effect on the explicitly described cluster atoms 
(especially the anions). 

In the Hartree-Fock calculations, / and g polarization 
functions as well as very diffuse functions had to be omit- 
ted for Ag and Au. The outermost spd exponents were 
reoptimized for the solid. (In the case of AuCl, two differ- 
ent basis sets were optimized for the two different crystal 
structures, respectively.) The CI basis set was taken from 
Ref. ^, again reoptimizing the diffuse exponents. The 
basis sets for the Hartree-Fock calculations are given in 
Table || 

III. HARTREE-FOCK CALCULATIONS 

As mentioned, HF self-consistent field (SCF) calcula- 
tions were performed with the program package crystal 
(Ref. ^). Results, together with those of the next section, 
are summarized in Table III. For AgCl, the SCF cohe- 
sive energy, Ecoh, is too small by 30% with respect to 
experiment. This, of course, reflects the well-known SCF 
error for bond-breaking. The situation is significantly 
improved for the lattice energy, Eiat, i-c. the energy for 
separating the crystal into Ag"*" and Cl~ ions — however, 
experiment is still underestimated by 20%. This has to 
be contrasted with the situation for the alkali halide crys- 
tal NaCl, which has a SCF lattice constant comparable 
to AgCl: the percentage of the SCF ermijs-for both Ecoh 
and Eiat are smaller by roughly 10-15%llJo. This should 
not be taken as an indication of increased covalent con- 
tributions in AgCl — the MuUiken population analysis 
shows a still nearly perfect ionic behaviour, with charges 
of ± 0.9 e for Ag and CI. However, it is certainly an 
indication that the correlation contribution of the Ag+ 
core is much more important than for the alkali halides. 
A similar conclusion can be drawn from considering the 
lattice constant: for AgCl, SCF yields a deviation from 
experiment of 0.5 A (9%) — in contrast to 'normal' 
quantum-chemical SCF results for covalent systems, it 
is an overestimation. This points to the importance of 
bond-shortening correlation effects like dynamic polar- 
ization and dispersion in AgCl. Of course, these effects 



are also acting in the alkali halides, but their magnitude 
is significantly smaller (3.5% in NaCl). Note that this is 
in line with the trend of the 'in-crystal' polarizabilities 
of the cations Na+ and Ag+ to be discussed below. The 
SCF bulk modulus of AgCl (~ 23 GPa) grossly underes- 
timates experiment, by more than 50%. 

Calculations for AuCl in the cubic rocksalt structure 
yield a SCF lattice constant and a SCF lattice energy 
both very similar to the corresponding values for AgCl. 
Also, according to our population analysis, the charges 
in AuCl are ±0.9 e, i.e. essentially the same as in AgCl. 
Thus, it may be concluded that the closed-shell repulsive 
potentials are quite similar for the (relativistic) Ag"*" and 
Au+ ions. The cohesive energy of AuCl, on the other 
hand, is smaller by more than 1 eV than in AgCl. This 
can be explained by the relativistic stabilization of the 
6s electron in the free Au atomE£l. Turning now to our 
results for the experimental tetragonal (lAi/amd) struc- 
ture of AuCl, we find it to be energetically nearly de- 
generate to the cubic structure, at the SCF level, with 
the latter very slightly lower (by ~1 mH). However, one 
can not conclude that at the Hartree-Fock limit (i.e. 
with a complete basis set) the experimental structure 
would still be higher in energy. In the rocksalt struc- 
ture, the d-occupancy on the chlorine atom is only 0.01 
e and omitting the d-function decreases the cohesive en- 
ergy by only 2 mH, whereas the d-population in the ex- 
perimental structure is 0.05 e and the energy decreases 
by 17 mH without this function. This results from the 
fact that a static d-polacization is symmetry-forbidden 
in the rocksalt structurellj. Thus, one might expect that 
additional polarization functions (e.g. / functions on 
Au) should further lower the energy of the experimental 
structure compared to the cubic one, and then already 
at the HF limit the experimentally observed structure 
would become more stable. As already indicated by the 
stronger static polarization effects, we find a tendency to- 
wards more covalent bonding in the experimental struc- 
ture: this shows up in the MuUiken population analysis, 
where the charge transfer Au — > CI is found to be not 
more than 0.5 e. In line with this change in bonding, 
the SCF overestimation of the lattice constant (5.8% for 
the experimental structure) is smaller than for the cu- 
bic one. (Note, at this point, that we performed only 
calculations for the isotropically expanded experimental 
structure, i.e. the ratio | was kept fixed as well as the 
parameter z — 0.19 which characterizes the z-component 
of the chlorine atoms with respect to the gold atoms.) 

For AgCl, Hartree-Fock andjiensity functional calcu- 
lations have been performedEHlEll. The Hartree-Fock re- 
sults are similar to ours (i.e. about 70 % of the cohesive 
energy is recovered and the lattice constant is overes- 
timated by 0.4-0.5 A). In Ref. ^ a detailed analysis 
of relativistic contributions was made, and reductions of 
the bond length, of about 1 %, and of the cohesive en- 
ergy, of about 15 mH, were found when effects of relativ- 
ity were included. Correlation-energy density-functional 
corrections improve the HF results significantly; however. 
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the lattice constant is still overestiniatedE3. An ab-initio 
treatment of electron correlation in silver and gold halides 
has not been attempted so far. 



IV. CORRELATION CALCULATIONS 

The calcu latia ns were done with the program pack- 
age MOLPROt^la. The correlated orbitals are 3s and 3p 
for CI; for Ag and Au we correlated the highest occu- 
pied d-shell (i.e. Ad for Ag, 5d for Au). We found in 
CCSD calculations for AgCl that additionally correlat- 
ing the inner-core 4s and 4p orbitals only slightly changes 
the results (e.g., the Ag-Cl increment for next neighbors 
changes from 0.008077 H to 0.008576 H). The incremen- 
tal expansion is made in terms of localized orbital groups 
which can be attributed either to the metal cations, M+, 
or are mainly centered on the chlorine anions. However, 
this does not imply any a priori assumption with regard 
to the ionicity of the system. 

In Tables and we give results for AgCl and AuCl 
(cubic structure) at two different lattice constants; cor- 
responding results for the tetragonal structure of AuCl 
are contained in Table VI. Let us first consider the 



intra-ionic 'one-body' contributions which are obtained 
when correlating only one embedded ion. As can be 
seen from the Tables, these energies are not too differ- 
ent from those of the free ions and therefore make only a 
minor contribution to the lattice energy. Much more im- 
portant are the van der Waals like inter-ionic two-body 
contributions; these are determined by correlating a pair 
of ions and subtracting off the corresponding intra-ionic 
one-body terms. We begin the discussion with results for 
the rocksalt structure. The dominant two-body contribu- 
tions clearly are the metal-chlorine terms. This is to be 
expected since these ions are nearest neighbors and the 
polarizability is not only high for the Cl~ ions but also 
for the metal ions due to their large number of electrons 
and large ionic radii. The chlorine-chlorine contributions 
are smaller by a factor of 4 ... 5 because of the larger 
distance between nearest-neighbor CI ions in the lattice, 
but in turn are still significantly more important (by a 
factor of 3 ... 4) than the metal-metal ones. It may be 
interesting to compare the AgCl correlation-energy incre- 
ments to those of NaCl which has about the same lattice 
constant. There, the Cl-Cl increments are the dominant 
ones, with the Na-Cl terms smaller by ~20 % and the 
Na-Na ones completely negligible — of course, this just 
reflects the relation between the metal-ion dipole polariz- 
abilities a: that of Ag+ being only a factor 2 smaller than 
for Cl~, Q!(Na"'"), on the other hand, being smaller than 
a(Cl^V-by more than an order of magnitude. It has been 
arguedcj that the difference in the strength of van der 
Waals interactions is responsible for the unusual features 
in the cohesive properties of AgCl as compared to the al- 
kali halides. In fact, according to our present CCSD(T) 
calculations (and those of Ref. Ilj), the ratio of intcr- 



ionic correlation contributions in AgCl and NaCl, at the 
equilibrium lattice constants, is '^4.4, which is not|-«ery 
different from the ratio of ~6 estimated by Buchero on 
the basis of a semiempirical fit to solid-state data. Com- 
paring now inter-ionic correlation in AgCl and AuCl, one 
notes an increase AgCl — > AuCl by about a factor 2. Of 
course, this has to do with the increase of the metal a 
but, as shown by the increase of the Cl-Cl increments, is 
also influenced by the reduction of the lattice constant 
in AuCl. In spite of the increased importance of inter- 
ionic van der Waals interaction in AuCl, the Au-Au in- 
crements are still rather small in absolute value: for a 
given Au-Au pair of nearest neighbors, the effect is only 
^0.04 eV. Let us now consider the experimental tetrag- 
onal AuCl structure. Here, distances between the ions 
are shorter, with a lower coordination number, however 
(smaller weight factors). This means that, for example, 
the Au-Cl two-body increment for next neighbors is much 
larger than in the rocksalt structure, but has to be mul- 
tiplied with a weight factor of two only (six in the cubic 
structure). The sum is still clearly dominated by the 
Au-Cl contributions which in total are roughly the same 
(within ^^15%) as in the cubic structure, the same holds 
for the Cl-Cl increments. In contrast to that, the corre- 
lation contribution between Au ions increases by a factor 
of '^3. The largest van der Waals attraction between an 
individual Au-Au pair (or, to be more precise, between 
the two c?^" shells) is now '^0.2 eV, i.e. relatively strong 
(comparable to the values given in Figure 36 of Ref. H 
at this distance of 3.2 A). However, this attraction alone 
would certainly not be strong enough to overcompensate 
the electrostatic repulsion of the positively charged ions 
(even if a reduced charge of 0.5 e, according to the Mul- 
liken population analysis, is assumed). We have no indi- 
cation, furthermore, for an additional stabilizing covalent 
interaction between the Au ions: the SCF overlap pop- 
ulation is -0.02. On the other hand, the 6s-like valence 
population on Au could give rise to non-negligible valence 
(and core-valence) correlation contributions — only, we 
cannot separate these contributions from those of the CI 
ions, since they cannot be attributed to different localized 
orbitals. 



In Table VII, we display 'in-crystal' ionization poten- 
tials (IP), dipole polarizabilities {a) and Cg coefficients 
for the cubic structures. The IP and a were calculated for 
an embedded single ion (treated with high-quality basis 
set and surrounded by pseudopotentials/point charges). 
For evaluating Cg , we multiplied individual 2-body incre- 
ments, for a given species, involving ions with the small- 
est internuclear distance in the lattice, by the sixth power 
of their distance. Of course, this implies the assumption 
of a pure van der Waals interaction between these ions, 
which is certainly only approximately satisfied. As al- 
ready mentioned, the in-crystal polarizability of Ag is 
larger by a factoij-ij^lO than that of Na, i.e. essentially 
comparable to Rbtj, and the a value of Au is even higher 
by a factor '^1.5. The Cg coefficients for Cl-Cl in AgCl, 
and even more so in AuCl, are larger than for the alkali 
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halidesB. While this trend can be quahtative ly ex plained 
for AgCl using the London formula, of. Table VII , the ra- 
tionalization for AuCl is more difficult. Apparently, the 
spill-over of charge on neighboring crystal ions is larger 
in AuCl, and the interaction is less van der Waals-like. 
The Cg coefficients for metal-chlorine interaction are next 
highest in magnitude. That of AgCl is in good agreement 
with RbCl, as expected from the similarity of the polar- 
izabilities, and the change from AgCl to AuCl roughly 
scales with the metal polarizabilities again. Turning now 
to the interactions between metal ions: our CCSD Cg 
value for Ag-Ag (71 a.u.) is in the range of 62 to 375 a.u. 
given in Ref. The Au-Au Ce coefficient (149 a.u.) is 
certainly very high compared to other metal-metal coef- 
ficients because of the high polarizability, but still signifi- 
cantly lower than both Au-Cl and Cl-Cl. This coefficient 
is also much lower when comparing to results from litera- 
ture for Cg coefficients characterizing the van der Waals 
interaction between two molecules, for example for the 
dimer (H3P-Au-Cl)2 a long-distance Cg limiting value of 
1830 a.u. has been calculatedQ. 

In order to check whether or not intra-ionic and 2-body 
inter-ionic correlation-energy increments give a reliable 
account of correlation effects in the systems considered, 
we also calculated several three-body increments. They 
were determined, as non-additivity corrections, in (em- 
bedded) clusters with three explicitly correlated ions, cf. 
footnotes to Tables 0, 0, and 0. As in earlier work 
they do not make an important contribution to the total 
energy. This does not seem to support the argument of 
Ref. ^ that three-body van der Waals interactions are 
important in AgCl. However, it is not clear whether the 
violation of the Cauchy relation of the elastic constants 
in AgCl would show up explicitly by means of 3-body 
increments in our formalism, or rather by an (angular) 
dependence of the 2-body ones on changes of lattice ge- 
ometry. 

Our final results are obtained by adding the (weighted) 
sum over all correlation-energy increments to the SCF en- 
ergies. By doing this for a range of lattice constants, we 
can monitor the influence of correlation on lattice con- 
stants and bulk moduli (see Table |l|). While HF cohe- 
sive energies underestimate the experimental values by 
30 % (AgCl) or 50 % (AuCl), the agreement is excel- 
lent at the correlated level for AgCl with deviations of 
less than 2%. For AuCl, we find a less good agreement 
with an overestimation of 34 mH (17 %) for the cohe- 
sive energy or 30 mH (7%) for the lattice energy. Lat- 
tice constants are strongly reduced due to the van der 
Waals interaction, by up to 0.8 A or 14 % (CCSD(T)) in 
the hypothetical cubic structure of AuCl. In the rocksalt 
structure, AuCl would have a shorter lattice constant (by 
0.3 A) than AgCl. The shortest Au-Cl distances found 
in the tetragonal structure of AuCl are 2.36 A, the short- 
est Au-Au and Cl-Cl distances are 3.22 A and 3.21 A , 
respectively (our calculated CCSD(T) values are 3.9 % 
shorter). The reduction of the distances compared to 
the cubic structure (CCSD(T) values: Au-Au = Cl-Cl 



= 3.69 A, Au-Cl = 2.61 A) may be explained with the 
rcduced-ipnic radii in systems with smaller coordination 
numbciO. The larger deviation of our results for the ex- 
perimental tetragonal AuCl structure as compared to the 
cubic AgCl one is probably due to the less good trans- 
ferability of the increments in the former structure (see 
footnotes to Tables and VI). When we extract the 

largest Au-Cl increment (which makes the most impor- 
tant correlation contribution to the lattice energy) from 
a cluster with three (instead of two) explicitly described 
ions, the cohesive energy is reduced by 12.5 mH and the 
deviation of the lattice constant from experiment shrinks 
to 3.5 %, at the CCSD(T) level. Note also that our cal- 
culations refer to zero temperature where AuCl is exper- 
imentally unstable. While the cubic structure for AuCl 
was slightly lower in energy (by 1 mH) than the experi- 
mental one, at the SCF level, the situation is reversed by 
correlation, and the experimental structure is calculated 
to be significantly more stable now (by 39 mH). About 
half of the stabilization is due to Au+-Cl^ pairs, the next 
important contribution (~40%) coming from the van der 
Waals interaction of the Au"*" cores. Both effects involve 
correlation of the Au 5c? shell. 



V. CONCLUSION 

We have extended the method of local increments to Ad 
and 5d transition- metal compounds. The results for co- 
hesive energy and geometry are in excellent agreement 
with experiment for AgCl and in good agreement for 
AuCl. For AuCl, we find the experimental structure 
with space group lAi/amd to be lower in energy by ~ 
39 mH than the hypothetical rocksalt structure. While 
the latter structure corresponds to nearly purely ionic 
bonding, the experimental structure features significant 
covalent bonding. The largest correlation contributions 
to the lattice energy come from pairs of nearest-neighbor 
metal/halide ions (up to >1 eV for a Au-Cl pair in the 
experimental structure of AuCl). However, while the 
metal-metal pairs give only very small van der Waals-like 
contributions in the cubic structure of AgCl and AuCl 
(<0.04 eV for a single M-M interaction), the correlation 
effect of Au-Au pairs reaches that of the Cl-Cl ones (~0.2 
eV per pair) in the tetragonal structure. The stabiliza- 
tion of the tetragonal structure with respect to the cubic 
one is essentially a correlation effect involving excitations 
from the Au 5d shells. The influence of correlation on 
the geometry turns out to be very important: the lattice 
constant is reduced by 9 % for AgCl and 10 % AuCl (ex- 
perimental structure), for AuCl in the rocksalt structure 
the reduction would be even larger (14 %). 



4 



ACKNOWLEDGMENTS 

KD and HS are grateful to Prof. P. Fulde (Dresden) 

for continous support. KD would like to thank Dares- 
bury Laboratory where this work was completed. 

* Present address: Daresbury Laboratory, Daresbury, 
Warrington, WA4 4AD, Great Britain 



^ Unkonventionelle Wechselwirkungen in der Chemie met- 
allischer Elernente: Bericht zum Schwerpunktprogramm 
"Neue Phanomene in der Chemie metallischer Ele- 
rnente mit abgeschlossenen inneren Elektronenzustanden" 
Deutsche Forschungsgemeinschaft, Editor: B. Krebs (VCH, 
Weinheim, 1991). 

2 H. Schmidbaur, Gold Bull. 23, 11 (1990). 

^ P. Pyykko, Chem. Rev. 97, 597 (1997). 

* P. Pyykko, N. Runeberg, and F. Mendizabal, Chem. Eur. J. 
3, 1451 (1997); P. Pyykko and F. Mendizabal, Chem. Eur. 
J. 3, 1458 (1997); P. Pyykko and F. Mendizabal, Inorg. 
Chem. (submitted). 

^ J. Li and P. Pyykko, Chem. Phys. Lett. 197, 586 (1992). 
^ E. M. W. Janssen, J. C. W. Folmer, and G. A. Wiegers, J. 

Less-Common Metals 38, 71 (1974); J. Strahle and K.-P. 

Lorcher, Z. Naturforsch. 29b, 266 (1974). 

E. M. W. Janssen, F. Pohlmann, and G. A. Wiegers, J. 

Less-Common Metals 45, 261 (1976). 

* R. Dovcsi, V. R. Saunders, and C. Roctti, M. Causa, 
N. M. Harrison, R. Orlando, and E. Apra, CRYSTAL 95 
User's Manual, Theoretical Chemistry Group, University 
of Torino (1996). 

® C. Pisani, R. Dovesi, and C. Roetti, Hartree-Fock Ah Initio 
Treatment of Crystalline Systems, edited by G. Bcrthier et 
al. Lecture Notes in Chemistry Vol. 48 (Springer, Berlin, 
1988). 

Quantum- Mechanical Ab-initio Calculation of the Proper- 
ties of Crystalline Materials, Editor: C. Pisani (Springer, 
Berhn, 1996). 

" H. Stoll, Phys. Rev. B 46, 6700 (1992); H. StoU, Chem. 
Phys. Lett. 191, 548 (1992). 

i^B. Paulus, P. Fulde, H. Stoh, Phys. Rev. B 51, 10572 
(1995); B. Paulus, P. Fulde, H. StoU, Phys. Rev. B 54, 
2556 (1996); S. Kalvoda, B. Paulus, P. Fulde, H. StoU, 
Phys. Rev. B 55, 4027 (1997); B. Paulus, F.-J. Shi, H. 
StoU, J. Phys.: Condens. Matter 9, 2745 (1997); M. Al- 
brecht, B. Paulus, and H. StoU, Phys. Rev. B 56, 7339 
(1997). 

" K. DoU, M. Dolg, P. Fulde, H. StoU, Phys. Rev. B 52, 4842 

(1995) ; K. DoU, M. Dolg, H. StoU, Phys. Rev. B 54, 13529 

(1996) ; K. DoU, M. Dolg, P. Fulde, H. StoU, Phys. Rev. 
B 55. 10282 (1997); K. DoU, M. Dolg, P. Fulde, H. StoU, 
Chemical Papers 51, 357 (1997). 

" K. DoU and H. StoU, Phys. Rev. B 56, 10121 (1997); Phys. 
Rev. B 57, 4327 (1998). 



D. Andrac, U. HauBermann, M. Dolg, H. StoU, and H. 
Preufi, Thcor. Chim. Acta 77, 123 (1990). 

D. E. Woon and T. H. Dunning, Jr., J. Chem. Phys. 98, 
1358 (1993). 

^'^ P. Fuentealba, H. Stoll, L. v. Szentpaly, P. Schwerdtfeger, 
and H. Preufi, J. Phys. B 16, L323 (1983). 
P. J. Hay W. R. Wadt, L. R. Kahn, and F. W. Bobrowicz, 
J. Chem. Phys. 69, 984 (1978). 

M. Prencipe, A. Zupan, R. Dovesi, E. Apra, and V. R. 
Saunders, Phys. Rev. B 51, 3391 (1995). 

E. Apra, E. Stefanovich, R. Dovesi, and C. Roetti, Chem. 
Phys. Lett. 186, 329 (1991); Chem. Phys. Lett. 197, 338 
(1992) (Erratum). 

21 N. J. M. Geipel and B. A. Hefi, Chem. Phys. Lett. 273, 62 
(1997). 

■^■^ MOLPRO is a package of ab-initio programs written by H.- 
J. Werner and P.J. Knowles, with contributions from J. 
Almlof, R. D. Amos, M. J. O. Deegan, F. Eckert, S. T. 
Elbert, C. Hampel, W. Meyer, A. Nicklass, K. Peterson, 
R. M. Pitzer, A. J. Stone, P. R. Taylor, M. E. Mura, P. 
Pulay, M. Schuetz, H. Stoll, T. Thorsteinsson, and D. L. 
Cooper. 

2^ H.-J. Werner and P.J. Knowles, J. Chem. Phys. 82, 5053 
(1985); P.J. Knowles and H.-J. Werner, Chem. Phys. Lett. 
115, 259 (1985); H.-J. Werner and P.J. Knowles, J. Chem. 
Phys. 89, 5803 (1988); P.J. Knowles and H.-J. Werner, 
Chem. Phys. Lett. 145, 514 (1988); C. Hampel, K. Peter- 
son, and H.-J. Werner, Chem. Phys. Lett. 190, 1 (1992); 
P.J. Knowles, C. Hampel, and H.-J. Werner, J. Chem. 
Phys. 99, 5219 (1993); M.J.O. Deegan and P.J. Knowles, 
Chem. Phys. Lett. 227, 321 (1994). 

2* M. Bucher, Phys. Rev. B 30, 947 (1984); Phys. Rev. B 35, 
6432 (1987). 

M. Wilson, P. A. Madden, and B. J. Costar-Cabral, J. Phys. 
Chem. 100, 1227 (1996). 
2® R. D. Shannon, Acta Cryst. A 32, 751 (1976). 

C. E. Moore, Atomic Energy Levels, NSRDS-NBS 35 / Vol. 
PHI, Nat. Bur. Standards (Washington, DC, 1949, 1952, 
1958). 

2* CRC Handbook of Chemistry and Physics, 75*'' edi- 
tion. Editor: David R. Lide (CRC Press, Boca Raton, 
1994/1995). 



5 



TABLE I. The first ionization potential (IP) for Ag and An is given in Hartrees units, evaluated at Hartree-Fock and 
correlated levels, in comparison with experimental values. 



Atom 


RHF 


CCSD 


CCSD(T) 


expt. 


27 




Ag 


0.233121 


0.269033 


0.272983 


0.27842 


Au 


0.282265 


0.326956 


0.331273 


0.33904 



TABLE IL Exponents and contraction coefficients (in parentheses) of the basis sets for the CRYSTAL calculation. In the case 
of AuCl, we optimized two different basis sets for the cubic and the experimental structure. 



Exponent (contraction) Exponent (contraction) 

AgCl AuCl 

|— I cubic lattice experim(M^tal lattice 

Ag Is, 2sp, 3spd pseudopotentiallfj Au Is, 2sp, 3spd, Aspdf pseudopotentiallfj 



4s 


9.0884420 (-1.9648132) 




5s 


20.1152990 (-0.1597614) 




7.5407310 (2.7332194) 






12.1934770 (0.7893559) 




2.7940050 (0.1991148) 






6.0396260 (-1.5714057) 


5s 


1.4801580 




6s 


1.3737210 


6s 


0.65 




7s 


0.65001 0.630 


7s 


0.16 




8s 


0.171 0.102 


4p 


4.4512400 (-6.0833780) 




5p 


8.6096650 (2.0982231) 




3.6752630 (6.4168543) 






7.3353260 (-3.0458670) 


5p 


1.2912880 (0.7539735) 




6p 


1.9129700 (0.3791452) 




0.6525780 (0.2730597) 






1.0576950 (0.6456428) 


6p 


0.38 




7p 


0.452 0.442 


Ad 


7.9947300 (-0.0163876) 




5d 


4.1439490 (-0.4058458) 




2.7847730 (0.2814107) 






3.5682570 (0.4275070) 




1.2097440 (0.4863264) 






1.3443240 (0.4755405) 




0.5053930 (0.3867258) 






0.5552890 (0.5610972) 


5d 


0.198851 




6rf 


0.192 0.188 


2sp, 3sp 


as in Ref. 


CI 


Is, 2sp, 3sp 


as in Ref. |l9| 


4sp 


0.308 




4sp 


0.314 0.312 


5sp 


0.113 




5sp 


0.115 0.108 


3d 






3d 


0.46 0.40 
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TABLE III. Final results for AgCl and AuCl. Energies are given in Hartree. 



AgCl 



AuCl (cubic structure) 



AuCl (exp. structure) 



HF CCSD CCSD(T) exp. 

E^oh 0.1421 0.1885 0.1997 0.20310 

Eic,t 0.2801 0.3315 0.3414 0-348fn 

a 6.00 A 5.57 A 5.52 A 5.51 AEl 

B 23 GPa 47 GPa 50 GPa 53.5 ± 0.5 GPi 



Ecoh 0.0994 0.1730 0.1936 

Eiat 0.2865 0.3738 0.3956 



a 



6.04 A 5.35 A 5.22 A 



aS 



Ecoh 0.0986 0.2067 0.2328 0.198SE!l 

Eiat 0.2858 0.4076 0.4348 0.4051 

a/c exp. +5.8 % exp. -3.1 % exp. -3.9 % 6.734 A/ 8.674 
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TABLE IV. Local correlation-energy increments (in Hartree) for cubic AgCl. The increments include weight factors appro- 
priate to one formula unit of the crystal. Position vectors (i,j,k) are given in units of o/2, where a is the lattice constant. 



lattice constant 

5.55 A 6. on A 





W'l >1 O it f;| cfi IT 


C'C'SD 


C'C'SD (T) 


(X'SD 


(.'(.' SI) (T) 


6iiib6d.tl6(d CI 


1 


-0.219821 


-0.228511 


-0 220706 


-0.229654 


free Cl~ 


.1 




-1-0 931 560 


-4-0 999073 


+0 931 560 


6iiibcclclccl Ag"'" 


1 


-0 344594 


-0 357267 


-0 344572 


-0 357224 




.1 


4-0 34.4551 


-1-0 3571 89 


-1-0 344551 


-1-0 357189 


C1(0,0,0)-C1(0,1,1) 


Q 


-0.010217 


-0.012170 


-0.006733 


-0.008127 




Q 

o 


-U.UUU400 


-U.UUU04D 


-u.uuuzyy 


-U.UUUoDz 


C1(0,0,0)-C1(2,1,1) 


12 


-0.000482 


-0.000580 


-0.000324 


-0.000393 


Ag(0,0,0)-Cl(l,0,0) 


6 


-0.048464 


-0.056076 


-0.032905 


-0.038437 


Ag(0,0,0)-Cl(l,l,l) 


8 


-0.001714 


-0.002040 


-0.001033 


-0.001234 


Ag(0,0,0)-Cl(2,l,0) 


24 


-0.000899 


-0.001069 


-0.000566 


-0.000675 


Ag(0,0,0)-Ag(0,l,l) 


6 


-0.002565 


-0.003000 


-0.001475 


-0.001724 


one-body-contribution to 




-1-0.002209 


-F0.002964 


-F0.001346 


-t-0.001864 


lattice energy 












two-body-contribution to 




-0.064769 


-0.075481 


-0.043335 


-0.050952 


lattice energy 












total correlation contribution 




-0.062587 


-0.072517 


-0.041989 


-0.049088 



to lattice energy" 



" We calculated two examples for three-body increments at 5.55 A: The increment C1(1,0,0)-C1(0,1,0)-C1(0,0,1) is 
-10.000041 H, at the CCSD level (CCSD(T): +0.000043): the CCSD increment for Ag(0,l,0)-Cl(0,0,0)-Cl(0,l,l) is 
-10.000103 H (CCSD(T): +0.000071 H). When taken from the latter cluster, CCSD results for one- and two-body 
increments are: Ag+: -0.344138 H, Cr : -0.219889 H, Ag-Cl: -0.048006 H, Cl-Cl: -0.011010 H; comparison with the 
corresponding values of the Table shows that the transferability of the increments is good. 
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TABLE V. Local correlation-energy increments (in Hartree) for cubic AuCl. The increments include weight factors appro- 
priate to one formula unit of the crystal. Position vectors (i,j,k) are given in units of a/2, where a is the lattice constant. 



lattice constant 

5.35 A 6. on A 





Wt'l^lit tiictcji' 


CX 'SD 


(:CSD(T) 


C'CSD 


('(:SD(T) 


6iiib6(dtl6<d CI 


1 


-0 216956 


-0.225227 


-0.218147 


-0.226611 


free CP 


.1 


+n 999073 


-1-0 931 560 


-1-0 999073 


4-0 931 560 


enibcddcd A.u"'" 


1 


-0 316662 


-0 329529 


-0 316643 


-0 329463 


11 CC iT. Li 


.1 


-i-f) 31 fifi98 


-1-0 399490 


-1-0 31 6698 


4-0 39Q490 


C1(0,0,0)-C1(0,1,1) 


Q 


-0.016820 


-0.020159 


-0.007916 


-0.009572 




Q 

o 


-U.UUiooy 


-U.UUiDU4 


-U.UUU4DU 


n f\f\f\KAQ 

-u.uuuo4y 


C1(0,0,0)-C1(2,1,1) 


12 


-0.000697 


-0.000849 


-0.000302 


-0.000365 


Au(0,0,0)-Cl(l,0,0) 


6 


-0.094949 


-0.110352 


-0.053367 


-0.062973 


Au(0,0,0)-Cl(l,l,l) 


8 


-0.003991 


-0.004798 


-0.001517 


-0.001818 


Au(0,0,0)-Cl(2,l,0) 


24 


-0.001851 


-0.002220 


-0.000767 


-0.000916 


Au(0,0,0)-Au(0,l,l) 


6 


-0.006699 


-0.007879 


-0.002897 


-0.003404 


one-body-contribution to 




-1-0.005083 


-F0.006224 


-F0.003911 


-t-0.004906 


lattice energy 












two-body-contribution to 




-0.126366 


-0.147861 


-0.067224 


-0.079597 


lattice energy 












total correlation contribution 




-0.121283 


-0.141637 


-0.063313 


-0.074691 



to lattice energy" 



" One three-body increment was calculated at 5.55 A: the CCSD contribution for Au(0,l,0)-Cl(0,0,0)-Cl(0,l,l) is 
-FO. 000233 H (CCSD(T): -F0.000172 H). CCSD results for one- and two-body increments, when taken from this 
cluster, are: Au+: -0.315803 (-0.316652) H, Cr: -0.217860 (-0.217040) H, Au-Cl: -0.075516 (-0.079530) H, Cl-Cl: 
-0.015054 (-0.013038 H). Comparison with the numbers in parentheses, which correspond to results from clusters with 
one or two explicitly described ions, shows again that the transferability is good. 
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TABLE VI. Local correlation-energy increments (in Hartree) for AuCl (experimental structure). Oexp represents the experi- 
mental lattice constants a=6.734 A and c=8.674 A. The increments include weight factors appropriate to one formula unit of 
the crystal. 



lattice constant 







0.95*a 


exp 


1.00*a 


exp 




weight factor 


CCSD 


CCSD(T) 


CCSD 


CCSD(T) 




1 




-0 228022 


-0 221 509 


-0 230664 


free Cl~ 


-1 


4-0 222073 


+0 231560 


-1-0 222073 


+0 231560 


PTnl^prlrlpr] An 


1 




-0 331 708 


-0 318760 


-0 332342 


frpp A n 


_! 


-1-0 316628 


-1-0 329420 


-1-0 316628 


-1-0 329420 


ClfO 25 IQVClfO 25 56) 


1/2 


-0 002949 


-0 003472 


-0 002257 


-0 002678 




X 


-0 004004 


-0 004789 


-0 00301 6 


-0 003697 


cun n 95 n iQVcun 5 n 95 -n nis'i 


2 


-0 006536 


-0 008034 


-0 004654 


-0 005809 


C1(0,0.25,0.19)-C1(0,0.75,-0.19) 


1 


-0.004097 


-0.004923 


-0.002514 


-0.003070 


01(0,0.25,0. 19)-C1(0. 5, 0.75, 0.06) 


2 


-0.003042 


-0.003814 


-0.002000 


-0.002556 


C1(0,0. 25, 0.19)-C1(0, 0.25,-0.44) 


1/2 


-0.000731 


-0.000950 


-0.000451 


-0.000596 


Au(0,0,0)-Cl(0, 0.25, 0.19) 


2 


-0.080882 


-0.092820 


-0.070456 


-0.081532 


AufO OVClfO 5 25 -0 06) 


4 


-0 022088 


-0 027312 


-0 016036 


-0 020052 


AiifTi (^^ Pirn n n AA^ 


9 




n 01 9799 




-\j.\J\JijO 1 \J 


Au(0,0,0)-Cl(0.5,0.25,0.31) 


4 


-0.002424 


-0.002952 


-0.001676 


-0.002044 


Au(0,0,0)-Cl(0,0.25,0.56) 


2 


-0.000482 


-0.000578 


-0.000342 


-0.000410 


Au(0,0,0)-Au(0.25,0.25,-0.25) 


2 


-0.014226 


-0.017044 


-0.009960 


-0.011908 


Au(0,0,0)-Au(0,0.5,0) 


1 


-0.004624 


-0.005527 


-0.003183 


-0.003799 


Cl(0,0.25,0.19)-Au(0,0.5,0)-Cl(0,0.75,-0.19)" 


1 


-h0.001329 


-h0.001231 


+0.001133 


+0.001075 


one-body-contribution to 




-h0.001255 


-h0.001250 


-0.001568 


-0.002026 


lattice energy 












two-body-contribution to 




-0.156385 


-0.184930 


-0.124153 


-0.147650 


lattice energy 












total correlation contribution 




-0.153801 


-0.182449 


-0.124588 


-0.148601 



to lattice energy 



" Following CCSD results for one- and two-body increments have been extracted from the above three-atom cluster at 
l.OOWeo^p: CI: -0.220213 H, Au: -0.316332 H, C1(0,0.25,0.19)-C1(0,0.75,-0.19): -0.003499 H, Au(0,0,0)-Cl (0,0.25,0.19): 
-0.061924 H. Comparison with the corresponding values of the Table shows that the transferability is less good than 
in the cubic structure, especially for the Au-Cl increment between next neighbors which makes the most important 
contribution to the lattice energy. 



TABLE VII. Ce coefficients, determined from CCSD two-body increments. In-crystal polarizabilities (a) and ionization 
potentials (IP) were calculated as in Ref. 14. The last column is a comparison with an estimate from London's formula. All 
quantities except for lattice constants are in atomic units. 



System lattice AE -Cs = AS x r« IPcat IPan a^at Oiau -f^^^t^ x AS 

constant a in A 

Ag-Cl (AgCl) 5^55 -0.008077 170 046 Oi 03 2L5 2^5 

Ag-Ag (AgCl) -0.000428 71 2.4 

Cl-Cl (AgCl) -0.001703 283 1.9 

Au-Cl (AuCl) 5.35 -0.015825 264 0.42 0.36 13.4 20.5 3.3 

Au-Au (AuCl) -0.001117 149 2.6 

Cl-Cl (AuCl) -0.002803 374 3.3 
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